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i— !■ Abstract 

O High energy inclusive hadron production in the central kinematical region is an- 

q^, alyzed within the models of unitarized pomeron. It is shown that the sum of multi- 

D | pomeron exchanges with intercept ap(0) > 1 reproduce qualitatively contribution of 

the triple pole (at t = 0) pomeron to inclusive cross section. Basing on this analogy 
we then suggest a general form of unitarized pomeron contributions (in particular the 
dipole or tripole pomeron) to inclusive cross section. They lead to a parabolic form 
qq ' of the rapidity distribution giving < n >oc In 3 s (tripole) or < n >oc In 2 s (dipole). 

0^ ■ The models considered with suggested parametrization of pj-dependence for cross sec- 



tions well describe the rapidity distributions data in pp and pp interactions at energy 
yfs > 200 GeV. The predictions for one particle inclusive production at LHC energies 
are given. 

o\ : 

Pomeron with intercept «p(0) = 1 + e, £ > is a very attractive model from the 
phenomenological point of view [1]. It gives a simple and compact parametrization for various 
high-energy soft processes (elastic and deep inelastic scattering, diffraction and others) and 
describes quite well a lot of data for high enough energy (for example, total cross sections 
and small-i (< 1 Gev 2 ) elastic scattering at y/s > 5 GeV). 

On the other hand at s — > oo the contribution of such a pomeron violates unitarity 
explicitly. The model leads to total cross section of hadron interaction behaving as cr t (s) oc 
(s/sq) 6 {sq = 1 GeV 2 ) in contradiction to the Heisenberg-Froissart-Martin-Lukaszuk theorem 

cr t <^\n 2 (s/s ). (1) 
m 2 

Thus the model can be considered only as phenomenological tool and must be improved in 
a some manner in order to restore unitarity, at least, to avoid a rough violation of unitarity 
bound on the total cross section. Let us remind that originally the contribution of of such a 
pomeron corresponds to a simple pole of partial amplitude in the plane of complex angular 
momentum. 

There are the several ways to restore unitarity. The most simple method to do that is to 
sum multipomeron diagrams up (Fig. [1]) 
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Figure 1: Multipomeron contributions to elastic scattering amplitude 



Starting on one pomeron exchange which can be written in the form 

a P (t) / \ a P (t) 



a(s, 



, t) = r) P {t)9ab{t) f-j = -g a (t)g b (t) f -i-j (2) 
where so = lGeV 2 and 

1 + exp(-iira P (t)) 

IP * - : 77\ ) 

— sirnrapyt) 

then going to the impact parameter representation 

oo 
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K S M = o — / dbbJ (gb)a(s,-q ), (3) 
8ns 

one can obtain under certing simplifying assumption (see below) the amplitudes 

H(s, b) = -1 J2 Ga{n)( * b{n) [-2ih(s, b)]-. (4) 

n=l 

It appears in the form Eq.(jl]) if we assume that two-hadrons-n-pomeron amplitude is propor- 
tional to the product of two-hadron-pomeron vertices (it is a pole approximation for inter- 
mediate states) as shown in Fig. [2J Moreover, assuming either G(n) = C n or G(n) = C n \fn). 
we obtain from Eq.(j4"]) two well known schemes of pomeron unitarization: eikonal [2] or 
quasi-eikonal [3] and quasi-£/-matrix models [4,5]. 



H(s,b) 



HC n Ci 



(1 - e~ 2lC ^ h ^) , if G ajb (n) = Cl b 



b 



(5) 

h(sA,) if G a , b (n) = C: b Vri. 



1 + 2iC a C b h(s, b) ' 



If ctp(t) = 1 + e + a'pt and g a ,b(t) = exp(5 a il t) one can find that at s — > oo in the both 
models 

af(s) w 8tteR 2 (s) ln(s/s ) w 8irea' P ln 2 (s/s ), (6) 

where R 2 (s) — B a + B b + a' P \n(s/so). 

This result provides a ground for another method for constructing the amplitude. The 
possible way is to consider partial amplitudes with more complicated singularities than usual 
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Figure 2: Amplitude of interaction of two hadrons with n pomerons in a pole approximation 
but with phenomenological factor G a {n). 



simple angular momentum poles just from the beginning. It is worth to emphasize that the 
factorization of residues is valid not only for the simple j-poles but also for any isolated 
j-singularity [6] of amplitude. Thus one can consider, for instance, double pole (dipole 
pomeron) [7] instead of simple pole. In this model cr t (s) oc ln(s/so) at s — > oo. Another 
interesting possibility comes from triple pole at t = (at t ^ because of unitarity it must 
be a pair of colliding at t — hard branch points producing a triple pole) [8] . Triple pomeron 
gives crt(s) oc ln 2 (s/s ). Both models lead not only to a very good description of the hadron 
total cross sections but the differential elastic cross sections, deep inelastic scattering and 
vector meson photoproduction as well. 

In the light of such a quite successful applications of unitarized pomeron models it is in- 
teresting to see how manyparticle processes can be described in this framework. We consider 
here one particle distribution in rapidity and pseudorapidity. In the next section we will re- 
mind some of the results obtained earlier for multipomeron exchanges. Then we suggest the 
specific parametrizations of p t and y dependence in the dipole and tripole pomeron models, 
compare them with the availble data for 200 < y/s < 1800 GeV and make predictions for 
future experiments at higher energy. 



1 Multipomeron exchanges in the model with op(0) > 1 

Due to generalized optic theorem differential cross section of one particle inclusive production 
(a + b — > c + X) in the central kinematical region where 

S = ip a + Pb) 2 -> OO, 

* = (Pa-Pc) 2 , \t\ -> OO, tU 2 2 2 

t \2 I —= m t =m c+Pl {') 

U = iPb-Pc) , > OO, S 

M 2 = (p a + Pb - Pc ) 2 , s/M 2 -> 1 

is related with the diagram of Fig. [3j 

More exactly, at large energy and for the simple pomeron pole with «p(0) — 1 = e 

^d 3 a ^ d 3 a ^ , ., , _ _. 

E—— = E- — — = 8irDisc M 2M(a + b + c -> a + b + c) = 
d A p dpid 2 p t e £ 

9a(0) l^) £ Vc(Pt) 

where E, pi,Pt are energy and momenta of the inclusive hadron c, g a ,b(to = 0) are the coupling 
vertices aPa,bPb, sq = lGeV 2 . 
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Figure 3: Pomeron contribution to inclusive production in the central region 



It is more convenient for what follows to use another set of variables, (pt,y) or (pt,rj), 

1 E I v 

y=-\n- — — , 77 = -ln(tantf/2), (9) 



where $ is the scattering angle of hadron c in the center mass system. With a rapidity 
variable the cross section Eq.([H]) is read as 



E 



d 3 



a 



d 3 p 



g a (0)e^-y) Vc ( p 2y(yo+y) gb ( 0) = g a (o) Vc ^ t )g b (o) e 2 ^. 



Here yo is the maximal value of rapidity y. In the c.m.s. 

- l/21n(s/m 2 ) <y< l/21n(s/m 2 ), i.e. y = l/21n(s/m t 2 ). 



(10) 



In the [9] the contribution to inclusive cross section of the diagrams given on the Fig. 
Hla have been calculated. It was shown that due to Abramovsky-Gribov-Kancheli rules [10] 
only input diagram with one pomeron exchange contributes, rest sum of diagrams vanishes. 
Thus the cross section is determined by Eq. ffTO"]) . 

Making use of the sum rule 



d 3 p d 3 a(ab — > cX) 
~E dfp 



=< n c > a t [s) 



(12) 



and taking into account that a t ~ a ln 2 (s/s ) at s — * oo one can find 

dn c 1 , w ~ , .„ (s/snY 



dy a t (s) 



ln 2 (s/s ) 



(13) 



where V c is the integral of v c (pf) over pf. Then integrating Eq. ffTBl over y we obtain a power 
growth for mean multiplicity, < n c >oc (s/so) £ / ln(s/so). Let us note that dn c /dy does not 
depend on y that is not supported by experimental data. 

However, more diagrams must be added to calculate the inclusive cross section under 
interest. These diagrams are shown on Fig. HJb an d were calculated in [11,12]. Likewise the 
case of diagrams Fig. HJa the sum of all contributions with the reggeons between hadrons a 
and b vanishes (see Appendix for details). 
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Figure 4: Multipomeron exchange dia- 
grams for one particle inclusive produc- 
tion, (a) - diagrams calculated in [9] , (b) 
- diagrams calculated in [11,12] 



Figure 5: The dominating contribu- 
tions to the central inclusive produc- 
tion at s — > oo 



Finally, due to the above mentioned cancelation at s — > oo the inclusive cross section in 
the central region is dominated by contribution of the diagrams of the Fig. [5] and can be 
written in a general form as 

E % = 9a ^^ yo ~ vKiPtWyo + y)g b (o). (n) 

If the input pomeron is simple j-pole and has intercept ap(0) = 1 + e then 

F(y ±y) = (y Q ±y) 2 . (15) 

It is necessary to note that this result exactly coincide with those which can be obtained if 
we assume from the beginning that pomeron at t = is the triple j-pole at j — 1. 

This fact and similar ones valid for elastic amplitude give us some ground for a more 
general assumptions and allows to consider the diagram of Fig. [3] with pomerons of arbitrary 
hardness at t — 0. 

If pomeron contribution to partial amplitude (of elastic scattering) at t = is proportional 
to l/(j — l) u+l [y < 2 because of the unitarity bound) then one can argue that in Eq. (|14p 

cLti I 

F(yo±y) = (yo±y) u and — oc — -(y -yY(y + y) u . (16) 

dy a t (s) 

We would like to remark that such a behaviour of dn/dy (at v > 0) is in a qualitative 
agreement with high energy experimental data, which show a rise dn/dy at y$ and a behaviour 
in y closed to parabolic form. Taking into account that such a pomeron leads to cr t (s) oc 
ln"(s/s ) ~ Vq one can find from Eq. ([T2"]) at s — ► oo. 

— (y = 0) oc \n u (s/s ) and < n > oc \n 1+u (s/s ). (17) 
dy 

It is known (see e.g. [13]) that a good description of the mean hadron multiplicity is 
achieved within a logarithmic energy dependence with v = 1 (dipole pomeron) or v = 
2 (tripole pomeron). The above mentioned properties of the unitarized pomeron models 
concerning of one particle inclusive distribution are rather attractive, therefore they should 
be checked out quantitatively with the data. We do that in the next section. 
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2 Comparison of the unitarized pomeron models with 
the data 



2.1 Experimental data 

Our aim is not the detailed description of all data, we would like to demonstrate only a 
possibility of the considered models to reproduce the main features of the high energy data. 
Evidently, at lower energy we need to add more Regge contributions which increase the 
number of the fitting parameters. To avoid an extra number of contributions and parameters 
we deal with the data on Ed 3 a/d 3 p at = 200, 540, 630, 900, 1800 GeV (240 points) [14,15] 
and on dn/drj [14,16] normalized to <Ji n (48 points). The data are shown on Fig.Q 
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Figure 6: Invariant cross section of 
charged hadrons pp — > charged + X, 
Ed 3 a / d 3 p. Data are taken from [14- 
16] 



Even for the chosen high energies we see 
a nontrivial dependence of cross sections 
on p t . One can clearly see that slope is 
changing with energy. Taking into account 
that pt-dependence in the pomeron contri- 
bution is coming only from vertex function 
v c(Pn), one can expect that the slope effect 
can be explained in the model only due 
to subasymptotic contributions. For ex- 
ample, in the dipole model it can be simple 
pole contribution with ap(0) = 1 as well as 
/-reggeon contributing at not highest en- 
ergies. Another important feature of the 
data is a changing of an exponential in- 
creasing Ed 3 a/d 3 p at small transverse mo- 
menta pt < 1 GeV for a power like increas- 
ing at higher p t (larger than 1 Gev). This 
transition can be parameterized in appro- 
priate way. The data on the Fig. [6] are 
known for < y >= while the models 
should be valid at an arbitrary y in the 
central region. 

Another set of the data, namely, dn/drj 
is more interesting and important for our 
aim. This observable can be obtained be 
integration of Ed 3 a/d 3 p over p t and with 
a transformation from y to rj. 



Pt,max Pt,7nax j 

dn 2ix f d 2 a 2ir f / m 2 d 2 a 

T = / A2 A P* d P* = / \ l 2 n 2 M A P* d P* 18 

di] u in J d 2 p t drj a in J W mf coslr y d 2 p t dy 

o o ' 

where <jj n is inelastic cross section (we fitted the data on dn/drj, presented in [14-16] just 
for for this case). In order to perform the integration one has to know the various vertex 
functions v c (p^) in d 2 a/d 2 p t dy which are not determined within any Regge model. Therefore 
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we parameterize them in a some form just to reproduce the experimental data. We would 
like to emphasize here that the explicit form of ^-dependence is not crucial for models under 
interest. It plays only a subsidiary role in obtaining da/dr] or dn/drj. Dependence of the 
differential cross section on y is more important for a verification of our approach. 



2.2 Double pomeron pole (or dipole) 

At y/s > 200 GeV we take into account dipole (d) and simple (p) poles with a(0) = 1 and 
/-reggeon for the both reggeons (upper and lower) in the diagram of Fig. [3j Thus, we have 
six terms for cross section 

d 3 a 

E dFp = 9ddVdd ( pt ^ y ° ~ y }( y ° + y) + 2 yo9d P Vdp{Pt) + 9ppV P p{pt) + 

9dfV df {pt)[{y - y)e £ i<y° + ^ + {y + y)e E f^-^}+ (19) 
9pfVpf(pt)[e e f {y °- y) + e £ f( yo+ ^] + 9ffVf f (p t )e 2e f ya 

where g a b, a,b = d,p, f are constants, £/ = «/(0) — 1, v a ^(pt) are vertex functions defined as 
follows 

v ab (Pt) = (1 + P 2 t /Plr» ab [z- Bpt + c(l + pVpIY"}- (20) 



2.3 Triple pomeron pole (or tripole) 

Generally there are four terms to be accounted for in this model. They are triple (t), 
double (d) and simple (p) pomeron poles together with /-reggeon. In total we have 12 terms 
including interferences. To avoid too many parameters we consider here a more simple model 
which has the form 

d 3 a 

E ~d?p = gttVtt ^ ZuZd + 9tpVt P (pt)(zu + z d ) + g P pVp P (pt) + 

g t fv tf (p t )[z u e E f(y°+y) + z d e £ f(y°-y)}+ (21) 

9pfVpf(Pt)[e £f{y °- y) + e £ f(y° + y^} + 9ffV ff (pt)e 2£ f yo 
where t- and c?-terms are combined to the single term as 

z u = P(yo - yf + (yo -y), z d = f3(y + yf + (y + y) 
and functions v a b(pt) are defined by Eq. f[2"Ul) . 



2.4 Simple pomeron pole 

This model (violating unitarity bound) we consider just to compare with the previous ones. 
We take into account two pomeron poles: one ("supercritical") has intercept «p(0) = 1 + e, 
another has a(0) = 1. 



d?a 

E ~d?p = 9ssVss ^ SuSd + 9spV S p{Pt){s u + s d ) + gppV pp {pt) + 

9sfV sf (pt) [s u e £ f^ +y ^ + S(te «7(w-i0] + ( 22 ) 
9pfVpf(Pt)[e Efiy °- y) + e E /^+«)] + g f fV ff (p t )e 2£ f yo 



with s u = e< y °-y\ s d = e £{yo+y) and v ab (p t ) taken from Eq.(|20l). 
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2.5 The data fit 



Parameters of the models as well as \ 2 obtained in the fits are given in the Table, description 
of the data is demonstrated on Figs. [7] and [HJ 



Table 1: Parameters of the models obtained from the data fit 



Parameters 


Dipole 


Tripole 


Simple pole 


9dd, gtti 9ss 


(mb) 1 / 2 


2.93 ± 0.006 


1.519±0.004 


817.9±0.2 


9 dp-, 9t%>-, 9sp 


(mb) 1 / 2 


-11.38 ± 0.018 


-8.061±0.019 


-1005.8±0.15 


9pp 


(mb) 1 / 2 


61.57±0.56 


79.54 ± 0.56 


1265.9±0.46 


9df,9tf,9sf 


(mb) 1 / 2 


0.61 ± 0.002 


1.002±0.002 


5.106±0.005 


9 P f 


(mb) 1 / 2 


-7 54±0 023) 


-28 39 ± 03 

Y7 • '.J ' ./ — 1 — V 7 • V 7 


-12.283±0.01 


9ff 


(mb) 1 / 2 


1175 9±7 2 

-1 LI . '.7 — 1 — 1 • — 


272 62 ± 50 


260 51±0 67 

^1 \J\J • \J X — 1 — \J • \J 1 


P 




031 ± 001 




e 






046 (< 001) 

\J « KJ^I V,/ 1 V7«V7V7X I 


e f 


305±0 005 


201 f<0 001) 


2±0 001 


Po 


GeV 


0.525±0.002 


0.518±0.002 


0.492±0.002 


Pi 


GeV 


0.070±0.001 


0.073±0.001) 


0.072±0.001 


B 


GeV" 1 


0.864±0.009 


0.931±0.008 


0.895±0.007 


A 


83.52±1.91 


82.85±1.80 


82.87±1.63 


n 


3.828±0.037 


3.872±0.037 


3.854±0.035 


^ddi fJ'ttj ^ss 


1.349±0.001 


1.352±0.001 


1.370 (< 0.001) 


H'dpj l^tpi f^sp 


1.289±0.001 


1.313±0.001 


1.378 (<0.001) 


l^pp 


2.066±0.012 


1.818±0.091 


1.417±0.001 


^df> fafi ^sf 


0.881±0.001 


0.388 (<0.001) 


0.257±0.002 


thf 


0.329±0.001 


0.469 (<0.001) 


0.314 (<0.001) 


Hi 


0.881±0.001 


0.597 (<0.001) 


0.713±0.001 


Xtot 


795.9 


809.4 


787.1 


X 2 /n.d.f. 


3.01 


2.98 


2.96 



One can see that theoretical curves in three models are very close each to other, at least 
for energies where data exist. It is not a surprise because the parameters of the tripole and 
simple models in fact mimic the dipole pomeron model. In the tripole model parameter f3 
is equal to 0.03 thus the terms containing (y — y Q ) 2 are not too important at the achieved 
energy, yfs < 1800 GeV. A similar situation occurs in the simple pomeron model where 
a strong cancelation among the s- and p-terms mimics a logarithmic behaviour at these 
energies. 

9ss(s/ s y + g pp w g ss + g pp + g ss e ln(s/ s ) 

However a difference between the models' predictions is increasing with energy. It can be 
seen clearly on the Figs. [9ll0|lll which demonstrate behaviour of dn(r] = 0)/dr], < n > and 
< p t > in energy. 

Concerning the average transferred momenta we would like to draw attention to the 
following property of the considered models. We have obtained lower values for < p t > 
than usually follows from a simple extrapolation of the data on Ed 3 a/d 3 p to p t = (see 
for example [15]). It is caused by the quite high values of Ed 3 a/d 3 p at p t — >• (Fig. [7]) 
which are necessary to obtain after integration over p t the correct values of dn/dr]. Thus the 
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Figure 7: prdependence of inclusive 
cross sections at high energies. Data 
are taken from [14,15]. Curves corre- 
spond to considered models. Red solid 
line - dipole pomeron model, blue long 
dashed line - tripole pomeron model, 
green doted line - simple pomeron pole 
with a(0) > 1. Predictions for three 
LHC energies are shown as well. 



Figure 8: Density of the produced hadrons 
as function of pseudorapidity and energy. 
Red (solid), blue (dashed) and green (dot- 
ted) lines are theoretical values corre- 
spondingly in dipole, tripole and simple 
pole pomeron model. Solid symbols cor- 
respond the data normalized to a in , open 
symbols correspond to data normalized to 
&nsd (not used in the fit procedure). Pre- 
dictions for three LHC energies are shown 
as well. 



considered models predict a fast increasing of the one-particle inclusive distribution at very 
small p t . 



3 Conclusion 

We have shown that the high energy experimental data on one-particle inclusive distribution 
can be described well in the models of unitarized pomeron, which do not violate unitarity 
Froissart bound for total cross section. The dipole (tripole) pomeron model, correspondingly 
lead to dn(y = 0)/dr] oc Ins (In 2 s) and < n >oc In 2 s (In 3 s). These models predict small 
differences in dn(y = 0)/drj and < n > at low LHC energies while they are increasing with 
energy. 
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Figure 9: Density of the produced 
hadrons at r] = as function of en- 
ergy. 
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Figure 10: Mean multiplicity of 
the produced hadrons as function 
of energy, calculated in the interval 
-3.5 < rj < 3.5. 
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Figure 11: Average transverse momentum of the produced hadrons as function of energy. 



4 Appendix 

Contribution to the discontinuity (AA4) in M 2 of the diagram with nix reggeons (/i of them 
are cut) in the upper part, with reggeons {I2 of them are cut) in the bottom part, with 
m reggeons (/ of them are cut) between hadrons a and b has the form 

m m\ ni2 
1=0 l 1= 0l 2 =0 

which can be rewritten as 

m mi m2 

8TiAM {m ' mi ' m2) (s } t } u) = £fl j2(-^) m ~ l (7)(- 1 ) mi ~ ll (T)(- 1 ) m2 ~ l2 {T) x 

1=0 h =0/2=0 1 

m _, mi _, m2 _, 

Y[2ImP(\t\,ki) EI 2ImP(| M |,A; n ) [1 2lmP(s,k i2 )G({k}) 

1=1 H=l «2 = 1 

(24) 
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where P(s,k) = i(s / 'so) ^ - *^ -1 is the reggeon propagator in (s, /^-representation, G({k}) 
contains all the vertex functions and the square brackets imply that the expression inside 
them must be integrated over the momenta of all reggeons provided the total momentum 
equals zero. 

Because of the relation 




m/0 



the Exp. (1241) leads to the following inclusive cross section 

rfS mi m.2 

e—= 87rA££-M (0 ' mi ' m2) ( M ,«) = £E j2 E(-i) mi -' 1 (T 1 1 )(- 1 ) m2 " /2 (? 2 2 ) x 

" P mi mi mi mi h=0 h=0 

mi mi 

I] 2ImP(|t|,A; n ) R 2lmP(\u\, k h )G a mi {k 1: k 2 , . . . , k mi )G b m2 (k[, k 2: . . . , k' m2 ) x 

.11=1 12 = 1 

(mi,mi),r t tl p „2\ 

(25) 

where p t is the transverse momentum of the inclusive particle c. 

Supposing that the vertex function v ^ mi ' m2 ^ either does not depend on ki, k\ or it has 
factored dependence on {ki} and {fc-} one can obtain for Ed 3 a/d 3 p an asymptotic behaviour 
similar to Eqs.dHD,(H5D. 
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